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Cuo.sCoo.sCb-FeCIs graphite intercalation compound is a three-dimensional short-range Ising 
spin glass with a spin freezing temperature T c (= 3.92 ±0.11 K). The stability of the spin glass 
phase in the presence of a magnetic field H is examined from (i) the spin freezing temperature 
Tf(uj,H) at which the differential field-cooled (FC) susceptibility 8Mfc(T, H) /dH coincides with 
the dispersion x'i^^T^H = 0) with the angular frequency cu, and (ii) the time dependence of zero- 
field cooled (ZFC) susceptibility xzfc after a ZFC aging protocol with a wait time tw (= 1.0 x 10 4 ). 
The relaxation rate S(t) (= dxzFc/dlnt) exhibits a local maximum at a characteristic time t cr , 
reflecting non-equilibrium aging dynamics. The peak time t cr (T,H) decreases with increasing H 
at the fixed temperature T (2.9 K < T < T c ). The spin freezing temperature Tf(cu,H) provides 
evidence for the instability of the spin glass phase in thermal equilibrium in a finite magnetic field. 
Both t cr (T,H) and Tf(u,H) exhibit certain scaling behavior predicted from the droplet picture, 
suggesting a dynamic crossover from SG dynamics to paramagnetic behavior in the presence of H. 

PACS numbers: 75.30.Kz, 75.40.Gb, 75.60.Lr, 75.50.Lk, 75.40.Mg 



I. INTRODUCTION 

It has been widely recognized that three-dimensional 
(3D) Ising spin glass (SG) with short range interactions 
undergoes a second order SG transition at least at zero 
magnetic field in thermal equilibrium at a finite spin 
freezing temperature T c } In spite of enormous amount 
of works, however, no consensus has been reached on 
the nature of the low temperature SG phase. Two dif- 
ferent pictures are proposed: mean-field (MF) picture 2 
and droplet pictured In the MF picture, the low tem- 
perature SG phase exhibits a replica symmetry break- 
ing. The SG ordered phase consists of infinite number of 
pure states organized in an ultrametric structure. In the 
droplet picture, in contrast, the SG phase is described 
by only two pure thermodynamic states related to each 
other by a global spin flip. 

The stability of the Ising SG phase in an external mag- 
netic field (H) has been one of the central issues. In the 
MF picture, the phase transition can survive in the pres- 
ence of low H, forming a critical line, so-called the de 
Almeida-Thouless (AT) line in the (ii", T) phase diagram: 

H(T) = A(l - T/T c y, (1) 

where A is a field amplitude and the exponent p = 3/2. 5 
This line separates the paramagnetic (PM) phase from 
the SG phase. There is no change of symmetry at the AT 
transition. The correlation length and relaxation time 
diverge on crossing this line. In the droplet picture, no 
phase transition occurs in the presence of even infinitesi- 
mal ii" as for a ferromagnet. So there is no AT line in the 
(ii", T) phase diagram. Any apparent transition would 
be an artifact related to the limited experimental time 
scale. Recent Monte Carlo (MC) simulations on the 3D 
Ising Edwards- Anderson (EA) 6 model with short range 
interactions suggest that the AT line is of dynamic origin 



and vanishes in the limit of infinite times. Takayama, 7 
and Takayama and Hukushima 8 have studied aging phe- 
nomena in the 3D Ising EA model under the ii'-shift ag- 
ing protocols, using nonequilibrium Monte Carlo simu- 
lations. They have shown that the characteristic length 
scales associated with the ii'-shift aging protocols exhibit 
unique scaling behavior. This scaling behavior implies 
the instability of the SG phase in the equilibrium limit 
even under an infinitesimal ii". Young and Katzgraber 9 
have also studied the possibility of the AT line for the 3D 
EA Ising SG using MC simulations. A finite-size scaling 
of the correlation length shows no indication of a transi- 
tion, in contrast to the zero-field case. This result also 
suggests that there is no AT line for the short range Ising 
SG. Experimentally Katori and ItoiS have reported the 
existence of the AT line with p = 1.49 in the 3D Ising SG 
Feo.sMno.sTiOs, supporting the MF picture. The criti- 
cal temperature T g (H) is determined as the one at which 
zero-field cooled (ZFC) susceptibility Xzfc starts to de- 
viate from the field-cooled (FC) susceptibility Xfc- The 
difference Ax = Xfc — Xzfc is a measure for the irre- 
versibility of the system. Mattson et al. 11 and Jonsson 
et al™ have shown from dynamic scaling analysis for the 
AC magnetic susceptibility of the same system that there 
is no SG transition in the presence of any ii", supporting 
the droplet picture. 

Cuo.sCoo.sC^-FeCls graphite bi-intercalation com- 
pound (GBIC) magnetically behaves like a 3D short- 
ranged Ising SG. 13,14 This compound undergoes a SG 
transition at T c = 3.92 ± 0.11 K in the absence of 
H. We have determined the spin freezing tempera- 
ture Tf{uj,H) at which the differential FC suscepti- 
bility dMpc(T^ H)/dH coincides with the dispersion 
x'ioo^T^H = 0), where oj (= 2nf) is the angular fre- 
quency. In the present paper we report experimental 
result on the time dependence of Xzfc after a ZFC 
aging protocol with a wait time tw (= 1.0 x 10 4 and 
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2.0 x 10 3 sec). We show that the relaxation rate S(t) (= 
dxzFc(^)/dln£) exhibits a local maximum at a charac- 
teristic time t cr , reflecting non-equilibrium aging dynam- 
ics. We find that t cr decreases with increasing H at the 
fixed T (2.9 K < T < T c ). The experimental results of 
Tf(uj,H) and t cr (T,H) are discussed in terms of predic- 
tions from both the MF picture and the droplet picture 
(see Sec. [HJ. We show that Tf(u,H) provides evidence 
for the instability of the SG phase in thermal equilib- 
rium in a finite H. We also show that both t cr (T,H) 
and Tf(u), H) exhibit certain scaling behavior, suggesting 
a dynamic crossover from SG dynamics to paramagnetic 
behavior in the presence of H. 

II. THEORETICAL BACKGROUND 
A. MF picture: equilibrium SG transition 

In the mean-field picture, each point on the AT line 
in the H-T phase diagram, T g (H), is a critical point 
of a continuous phase transition in the presence of 
H which exhibits critical divergence of the correlation 
length £ C (T, H) and the correlation time r c (T, H) on ap- 
proaching T g (H) from the high-temperature side. 

T C (T,H) = TS(UT,H)/L y, (2) 

or 

t c (T,H)=t5[T/T 9 {H)-1]-*, (3) 

where 

UT,H) = L [T/T g (H)-l}-», (4) 

To and Lq are microscopic units of time and length, and 
z and v are critical exponents (x = vz). From the onset 
of the out-of-equilibrium behavior of the AC susceptibil- 
ity, the freezing temperature Tf(w,H) can be extracted, 
where the probing time t is given by t = 2tt/uj for the 
AC susceptibility. Then the critical temperature T g (H) 
is derived from the condition r c (Xy, H) = t = 2tt/lj^ 

T g (H) = T^H)/^/^)- 1 ^ + 1]. (5) 

If the AT line exists in the system, it is required that 
T g (H) can be uniquely determined using the appropriate 
values of Tq and x. Note that Tq and x are assumed to 
be the microscopic time and the dynamic critical expo- 
nent at H = 0, respectively. The line T g (H), which is 
independent of is assumed to form a AT line of the 
second order in the H-T plane, where T g (H = 0) = T c . 

B. Droplet model: scaling form for aging dynamics 

The situation is rather different in the the droplet 
model. After the SG system is cooled to T (< T c ) 



through the ZFC aging protocol, the size of domain de- 
fined by Rr(ta) grows with the age t a from t a = and 
reaches Rr(tw) just before the field H is turned on at 
the observation time t = t a — tw = (or t a = tw)- 8 Here 
Rr(t a ) is expressed by 

RT{ta) = L T {t a /T* ) bT l T % (6) 

where Lt is a characteristic length, and b is an exponent. 
The aging behavior in Xzfc is observed as a function of t. 
After t = 0, a probing length Rt r (t) corresponding to the 
maximum size of excitation grows with t, in a similar way 
as Rrita)- Here Rt r (t) is the mean size of the subdomain 
in the transient regime and is expressed as£ 

R tr {t) = L T {t/T* ) bT l T c. (7) 

The quasi equilibrium relaxation occurs first through 
local spin arrangements in length Rt r (t) <C Rr{tw)i 
followed by non-equilibrium relaxation due to domain 
growth, when Rt r (t) ~ Rritw), so that a crossover be- 
tween the short-time quasi-equilibrium relaxation and 
the non-equilibrium relaxation at longer observation 
times is expected to occur near t « tw I n the time 
range < t < t cr , which we call the transient regime of 
i^-shift processes, SG domains in local equilibrium with 
respect to H > grow within the SG domains which were 
in local equilibrium with respect to H = at t a = tw 

Here we introduce two characteristic lengths Rw and 
R cr defined by 

Rw = Rr(t a = tw) = L T (t w /rS) bT/T % (8a) 

Rcr = Rtr(t = t cr ) = L T (t cr / T*f T ^ , (8b) 

and the magnetic field crossover length Rh defined by 12 

RH^lnKl-T/T^ffH- 1 } 6 , (8c) 

where lu is a constant associated with a characteristic 
length, S = (d/2 - 6)~ l and a eff = Ov - (3/2. The defi- 
nition of <i, /?, and v is given in the previous paper. 13 
According to Takayama and Hukushima, 8 it is predicted 
that Y (= Rcv/Rh) is described by a scaling function of 
X(=R W /R H ): 

Y = R cr / Rh = f(Rw/R H ) = f(X). (9) 

The H-shih aging process is nothing but a dynamic 
crossover from the SG phase in H = to the param- 
agnetic state in H > in the equilibrium limit. In the 
limit of large Rw/Rh, Rcv/Rh becomes constant. 

For the AC magnetic susceptibility measurement with 
the angular frequency a;, the size Rt r (t = 2tt/co) is 
the mean size of droplets that can respond to the AC 
field of frequency uo. Then the onset of the SG out-of- 
equilibrium dynamics, T/(t, H), can be related to Rh by 
the condition^ 

Rh = Rtr(t = 2tt/cj), (10a) 

or 

(l-T f (w,H)/T c ) a "f/H = {LT/lnfl^n/w^f^l^fl 5 . 

(10b) 
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III. EXPERIMENTAL PROCEDURE 

The detail of sample characterization and sample 
preparation of Cuo.sCoo.sC^-FeCla GBIC was provided 
in our previous papers^AA The DC magnetization was 
measured using a SQUID magnetometer (Quantum De- 
sign, MPMS XL-5) with an ultra low field capability op- 
tion. The remnant magnetic field was reduced to zero 
field (exactly less than 3 mOe) at 298 K. The time (t) 
dependence of the zero-field cooled (ZFC) magnetiza- 
tion (Mzfc) was measured at fixed points in the (H,T) 
plane. The following ZFC aging protocol was made be- 
fore the measurement. First the sample was annealed at 
50 K for 1.2 x 10 3 sec in the absence of H. Then the 
system was quenched from 50 K to T (< T c ). It was kept 
at T for a wait time tw (= 2.0 x 10 3 sec and 1.0 x 10 4 
sec). After the wait time, an external magnetic field H 
is applied along any direction perpendicular to the c-axis 
at t = (or t a = tw): where t a is an age and t is an 
observation time (t a = t + tw)- The measurements of 
Xzfc vs t were carried out at the fixed T (2.9 < T < T c ) 
when H is varied as a parameter: 1 < H < 400 Oe. 

Under the above conditions, XzFc(t) (= MzFc(t) / H) 
was measured as a function of t. The relaxation rate 
defined by S(t) (= dxzFc(t) /dint) exhibits a peak 
at a characteristic time t cr . 15 Theoretically 16,17 and 
experimentall}ii2ii2i2£l it has been accepted that the time 
variation of XzFc(t) may be described by a stretched 
exponential relaxation form: 

XzFc(t) =Xo~ B t~ m expHi/T) 1 - 11 ], (11) 

where the exponent m may be positive and is very close to 
zero, n is between and 1, r is a characteristic relaxation 
time, xo, ^? and Bo are constants, and n is a stretched 
exponential exponent. The least-squares fit of the data of 
Xzfc vs t to Eq.(ll) with m = yields the parameters 
n, r, and Bq as a function of H and T in the H-T plane. 



IV. RESULT 

A. Line T g o(H) from the measurement of dAx/dT 
vs T 

In our previous paper— we have discussed the H de- 
pendence of the spin freezing temperature in detail. The 
result is summarized as follows. Figure da) shows the 
T dependence of the difference between Xfc and Xzfc 
(Ax = Xfc — Xzfc) at various H in the vicinity of T c . 
The transition temperature is defined as a temperature at 
which Ax = 0. Since it is difficult to determine exactly 
the transition temperature from Fig. ^a), we redefine 
the line T g o(H) as one at which dAx/dT vs T shows a 
local minimum for each H (see Fig. [Jb)). In Fig. [21 we 
show the H dependence of T g o(H) thus obtained. The 
least squares fit of the data of H vs T (or H vs T g o(H) 
in Fig. [21 to Eq.fTJ yields the parameters p and A, where 
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FIG. 1: (Color online) T dependence of (a) A% = xfc — Xzfc 
and (b) dAx/dT at various H for Cuo.sCoo.sCb-FeCb GBIC. 



T c is fixed as T c = 3.92 ±0.11 K. The exponent p is 
equal to 1.52 ± 0.10, which is close to the AT exponent 
= 3/2). The value of A is A exp = 3.38 ±0.67 kOe. For 
convenience, hereafter this critical line is denoted as the 
line T g o(H) (dAx/dT vs T). In the classical q vector SG 
model, the value of A is predicted as 18 - 21 

Atheory = [S/(q ± l)(q ± 2)] x ' 2 (k B T c /gfi B S) , (12) 

where S is the spin, g is the Lande ^-factor, and q is the 
number of spin degrees of freedom (q = 1 for Ising SG and 
g = 3 for Heisenberg SG). Here we assume g = 2. The 
value of £ is estimated as 5 = 1.56 for our system with 
the stoichiometry C 5 .26(Cuo.5Coo.5Cl2)o.47(FeCl 3 )o.53, 13 
where 5 = 1/2 for Co 2+ and Cu 2+ and S = 5/2 for Fe 3+ . 
Then the value of A t heory is estimated as 21.6 kOe for 
q = l (Ising symmetry), leading to A exp /A theory = 0.16. 
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FIG. 2: (Color online) H-T phase diagram, (i) The local 
minimum temperature of dAx/dT vs T at H (♦) [denoted 
as the line T g0 (H) (dAx/dT vs T)]. A solid line is a least- 
squares fitting curve which is expressed by Eq.Q with T c , p 
and A given in the text, (ii) Tf(cj, H) vs iif between / = 0.01 
Hz and 1 kHz. Tf(u, H) is defined as a temperature at which 
Xfc(T, H) crosses x'(a;, T, if = 0) at each cj. 



These results are similar to those reported by Katori and 
Itc 10 for Feo.5Mno. 5 Ti0 3 : p = 1.49 and A exp /A theory = 
0.52 (A exp = 42.6 kOe, A theory = 82.1 kOe). Although 
such result apparently supports the AT theory, it does 
not provide any evidence for an equilibrium SG transition 
in the presence of H. 



B. Line Tf(u,H) where 
dM FC (T, H)/dH = xV, T, H = 0) 

In our previous paper— we have discussed the 
spin freezing temperature Tf(u,H) of Cuo.5Coo.5Cl2- 
FeCl3, at which dMpc{T, H)/dH crosses the dispersion 
X r {w,T,H = 0) at each /. Such definition of Tf(uj,H) 
has been also used by Mattsson et al. 11 and Jonsson 
et al. 12 for T f (u,H) in Fe .5Mn .5TiO 3 . In Fig. we 
show the plot of Tf(uj,H) as a function of H where / 
(0.01 < / < 1000 Hz) is varied as a parameter. The line 
Tf(uu,H) is strongly dependent on /. If Tf(u,H) In the 
limit of / = coincides with the line T g o(H) (dAx/dT 
vs T) as shown in Fig. [2J then this line would be an 
equilibrium AT line. 

It has been reported that the SG phase transition oc- 
curs only at zero field in Feo.sMno.sTiOsiA^ The disap- 
perance of the AT line at finite fields is consistent with 
the prediction derived from the droplet picture. By using 
the same method (see Sec.[Hj), here we examine whether 
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FIG. 3: (Color online) Plot of T g (H) (= T f (u,H)/[l + 
(2tv /ujtS)~ 1/x ]) vs H with x = 10.3 and r * = 5.27 x 10" 6 
sec, where cu (= 27r/) is the angular frequency. 



the system exhibits the critical slowing down in the pres- 
ence of H. If the SG transition occurs at T g (H) in ther- 
mal equilibrium as is predicted from the MF picture, it 
is required that the critical temperature T g (H) evaluated 
from Eq.(0) should fall on a unique curve T = T g (H), 
which is independent of lj. 

We calculate the line T g (H) using Eq.© with the data 
of Tf(u,H) for each lj. Here we use Tq = 5.27 x 10" 6 
sec, T c = 3.92 K, and x = zv = 10.3, which are obtained 
from the analysis of the absorption ^'(o^T, = 0) in 
the previous paper i& Figure 01 shows the line T g (H) thus 
obtained. We find that the line T g (H) is still strongly 
dependent on lj, although the line T g (H) tends to ap- 
proach the line T g o(H) (dAx/dT). We note that the line 
T g (H) is still dependent of lo even if x is equal to 20, 
which is unphysically large. This result suggests that no 
equilibrium SG phase transition occurs in the presence 
of H, which is inconsistent with the prediction from the 
MF picture. 

Now we examine the validity of the droplet picture. 
As is described in Sec. |H] the onset of the SG out- 
of-equilibrium dynamics, Tf(t,H), is related to Rh by 
Eq.flBB with X 1 = (t/rS) T f^ H)/Tc and t = 2h/ljM 
When Yi is defined as Y x = (1 - T f (t, H)/T c ) a ^f / H , 
Eq. ([10b|) is rewritten as 



y 1 « e^xi' 5 , 



(13) 



with ( = Lt/Ih- Figure El shows the plot of ln(Yi) as a 
function of ln(Xi) with the best choice of a e ff = 0.55 ± 
0.05, where Tq = 5.27 x 10 -6 sec. When ony the data with 
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FIG. 4: (Color online) Scaling plot of ln(Yi) vs ln(Xi) 
for the data with T f (u,H)/T c < 0.9, where Yi = (1 - 
T f {uj,H)/T c ) a -ff /H and Xi = (27r/cjr * ) T ^ (t '^ )/Tc . The best 
result for the scaling plot is obtained when a e // = 0.55, where 
T c (= 3.92 K) and r * (= 5.27 x 10~ 6 sec) are fixed. The least- 
sqaures fitting curve (Yi ~ / ' Y ' J/d 
line with a slope of b/5 = 



^ C^X^ ) is denoted by a straight 
0.194 and C = 1.5 X 10" 4 (5 = 0.81). 



Tf(t,H)/T c < 0.9 are used, almost all the data fall well 
on an unique straight line with the slope b/5 = 0.198 ± 
0.002 and ln(C 1/(5 ) = -10.85 ± 0.02. Using the value of 
b = 0.16 for x"(a;,T = 3.75 K,t) with / = 0.05 Rz^ 5 
can be estimated as 5 = 0.81. This value of 5 is slightly 
deviated from the calculated value of 5 using the relation 
5 = (d/2 - 0)- 1 = 0.73 with d = 3 and = 0.13. Our 
value of a e f f is much larger than the value predicted from 
the relation (a eff = 0V-/3/2 w 0), where (3 = 0.36±0.03, 
v = 1.4 ± 0.2, and = 0.13 ± 0.02 which are obtained in 
the previous paper. 13 Using 5 = 0.81, the value of £ is 
calculated as (1.5±0.1)xl0 -4 . This value will be used in 
Sec. It is interesting to compare our result with that 
reported by Jonsson et al. for Feo.sMno.sTiOs. 12 They 
have reported the values of a e ff — 0.45 and b/5 = 0.17 
for the high-temperature (HT) data with Tf(cu,H)/T c < 
0.95 and H < 5 kOe, and a eff = 0.25, b/5 = 0.143 for the 
low-temperature (LT) data, where 6 = 0.2 and 5 = 0.77. 
Our result (a e // = 0.55 and b/5 = 0.198) is close that 
obtained from the HT data. 



C. Stretched exponetial relaxation form xzFc{t) at 
H = 1 Oe 



We have measured the t dependence of XzFc(t) af- 
ter the ZFC aging protocol described in Sec. IIIII Fig- 
ures Efa) and (b) show the t dependence of XzFc(f) 
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FIG. 5: (Color online) t dependence of (a) AxzFc(t) [= 
XzFc(t) — xzFcit — 0)] and (b) the relaxation rate S(t) 
(= dxzFc(t)/d\nt) at various T (3.0 < T < 4.7 K). t w = 
2.0 x 10 3 sec. if = 1 Oe. 



and the relaxation rate S(t) at various T, respectively, 
where tw = 2.0 x 10 3 sec and H = 1 Oe. The relax- 
ation rate 5(t) exhibits a peak at t = t cr , which shifts to 
the short-t side with increasing T. The ZFC susceptibil- 
ity XzFc(t) is well described by Eq.JHJ with m = for 
1.0 x 10 2 < t < 6.0 x 10 4 sec. The least-squares fit of the 
data of Xzfc vs t yields the parameters n and r for each 
T. In Fig.UJa), we show the T dependence of t cr and r. 
Both £ cr and r increase with decreasing T: r = 3.4 x 10 4 
sec at T — 3 K. Note that r is larger than t cr below T c . 
In the inset of Fig. a) we show the plot of l/t cr and 1/r 
as a function of 1/T. The least-squares fit of the data of 
1/r vs 1/T for T < T c to the formic 



1/r = ci exp(-c 2 T c /T), 



(14) 



yields the parameters c\ = 19. 7 ± 11.5 sec 1 and c 2 = 
10.82 ± 0.56. For the data of l/t cr vs 1/T, similarly we 
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FIG. 6: (Color online) (a) T dependence of t cr and r. t cr is 
a characteristic time at which S(t) exhibits a peak and r is 
the relaxation time for the stretched exponential relaxation. 
H — 1 Oe. tw — 2.0 x 10 3 sec. The inset shows the plots 
of I /t cr vs 1/T and 1/r vs 1/T for T < T c . The solid lines 
are least-squares fits to Ea. (|14f> for both l/t cr vs 1/T and 
1/r vs 1/T. The fitting parameters are given in the text, (b) 
T dependence of the stretched exponential exponent n. The 
solid line is a guide to the eyes. 
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FIG. 7: (Color online) t dependence of (a) xzFc(t) and (b) 
S(t) at various H (20 < H < 300 Oe). t w = 1.0 x 10 4 sec. 
T 3.50 K. 



obtain c x = 1.01 ± 0.48 sec -1 and c 2 = 7.81 ± 0.46. Note 
that our value of C2 is much larger than those derived 
by Hoogerbeets et al. 18 (c2 = 2.5) from their analysis 
of thermoremnant magnetization (TRM) relaxation mea- 
surements on four canonical SG systems: Ag: Mn (2.6 
at. %), Ag: Mn (4.1 at. %), Ag: [Mn (2.6 at. %) + Sb 
(0.46 at. %)], and Cu: Mn (4.0 at. %). Note that in their 
work the stretched exponential is taken as representative 
of the short time (t < tw) relaxation. In Fig. EJb) we 
show the T dependence of n. The exponent n exhibits a 
local maximum around T = 3.3 K (T/T c = 0.84), a local 



minimum around 3.5 - 3.7 K (T/T c = 0.89 - 0.94), and 
a sharp peak just above T c . It drastically decreases with 
further increasing T above T c . The overall variation of n 
below T c is similar to that reported by BontempsSi and 
Bontemps and Orbach 23 for Eu .4Sr . 6 S (T c = 1.50 K): 
n has a local maximum at T/T c = 0.87, a local minimum 
at T/T c = 0.91, and a local maximum at T = T c . In 
summary, the transition between the SG and PM phases 
at H = 1 Oe is characterized by the peak of n and the 
drastic increase of t cr (or r) with decreasing T around 
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FIG. 8: (Color online) t dependence of (a) xzFc(t) and (b) 
S(t) at various H (150 < H < 400 Oe). t w = 1.0 x 10 4 sec. 
T = 3.10 K. 



D. if and T dependence of £ cr and r 

We have also measured the t dependence of xzFc(t) 
after the same ZFC aging protocol described in Sec. Hill 
where t w = 1.0 x 10 4 sec, T = 2.9, 3.0, 3.1, 3.3, 3.5, 
3.65, and 3.75 K and H = 20, 50, 100, 150, 200, 300, 
and 400 Oe. Figures Efa) and (b) show the t dependence 
of XzFc(t) an d S(t) at T = 3.5 K, respectively. Figures 
EI a) and (b) show the t dependence of xzFc(t) an d S(t) 
at T = 3.1 K, respectively. Here H is chosen such that 
the point at (T,H) crosses the line T g o(H) (dAx/dT vs 
T) in the H-T plane. The relaxation S{t) has a peak at 
t = t cr , which shifts to the low-t side with increasing H. 
The peak height of S(t) at £ = t cr drastically decreases 
with increasing H. The least-squares fit of the data of 
XzFc(t) vs t for 1.0 x 10 2 sec < t < t w to Eq.fTTJ with 
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FIG. 9: (Color online) Plot of t cr vs r for various H between 
50 and 400 Oe at each T. A solid line is a least-squares fitting 
curve denoted by x cr — t cr /r = 0.80. tw = 1-0 X 10 4 sec. 



m = yields the exponent n and r at each H and T, 
where r is the relaxation time of the stretched exponen- 
tial relaxation. Note that the data (xzfc vs t) are not 
fitted well with Eq.JHJ with m ^ 0: the value of n is 
dependent on the small change in m. 

The values of r and t cr have been determined from 
the t dependence of xzFc(t) an d S(t) at fixed H and T. 
In Fig. H we make a plot of t cr vs r at each T (2.9 < 
T < 3.75 K) for measured values of H. We find that t cr 
is proportional to r in the range < r < 6.0 x 10 3 sec: 
x cr = t cr /r = 0.80 ±0.01. This ratio is different from the 
ratio (x cr = 1) predicted from Ea.([TT]) with m = 0, which 
is independent of n. A slight deviation from this linear 
relation is observed for r > 6.0 x 10 3 corresponding to the 
case of low H (H < 50 Oe), partly because of the larger 
uncertainty in r. Since the H and T dependence of t cr 
is essentially the same as that of r, hereafter we discuss 
only the H and T dependence of t cr . Figure EJa) shows 
the H dependence of t cr at fixed T (2.9 <T< 3.75 K) 
below T c . At higher T, the decrease of t cr with H occurs 
at lower H. Figure ITUT b) shows the T dependence of t cr 
at fixed H below T c . 



V. DISCUSSION 

We examine the validity of the scaling relation given by 
Eq.Q using the data of t cr (T,H) with t w = 1.0 x 10 4 
sec in Figs. ITUT a) and (b). Note that t cr (T,H) is the 
peak time of S(t) vs t for H 2.9 < T < 3.75 K and 
5 < H < 500 Oe. We define the following two parameters 
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FIG. 10: (Color online) (a) t cr vs H at various T and (b) t cr 
vs T at various H. The solid lines are guides to the eyes. 

X and Y given by 

X = C[H 5/b {t w /r^) T/Tc (l - T/T c )-(W a ^f}\ (15a) 
and 

Y = C[H 5/h (t cr (T, #)/r*) T/Tc (l - T/T c )-^ s ^ff}\ 

(15b) 

respectively. In Figs. lTTT a) and (b) we show a scaling plot 
of Y vs X with r * = 5.27 x 1(T 6 sec and T c = 3.92 K. If 
we use the parameters S/b (= 5.05), a e ff (= 0.55), and £ 
(= 1.5 x 10 -4 ), which are determined independently from 
Fig. 01 we find that almost all the points (X, Y) fall well 
on a single curve. Here we use the value of b (= 0.16). 13 
In fact, the scaling form of Figs. ITTT a) and (b) is not so 
sensitive to the choice of b for 0.14 < b < 0.18. Similar 
scaling curve has been reported by Jonsson et al^2& 
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FIG. 11: (Color online) (a) Scaling plot of 7 vs I. X and 
Y are denned by Eas. (|15af> and (|15b|> . respectively, where 
S/b = 5.15, b = 0.16, a e ff = 0.55, ( = 1.5 x 10" 4 , T c = 3.92 
K, and Tq — 5.27 x 10 -6 sec. t cr is a time of the peak for 
S(t) vs t with t w = 1.0 X 10 4 sec for 2.9 < T < 3.75 K and 
5 < H < 500 Oe. (b) The same scaling plot of Y vs X for 
< X < 0.6. The solid lines are denoted by Ea. (fT6j) with c 
= 0, 0.15, 0.3, and 0.71 where 5 = 0.81. Note that the solid 
line with c = 0.71 also coincides with that derived from the 
least-squares fit of the data to Eq.flSJ with c = 0.71 =b 0.06 
and 8 = 0.81 ± 0.07. 



from the aging behavior of xzFcif) f° r Ag 11 at. % Mn 
under the T-shift protocol. 

Our curve is well described by a straight line given by 
Y = X for X < 0.1. It deviates from the straight line 
with increasing X and tends to saturate around X w 1.6. 
There is a gradual crossover between accumulative and 
chaos around X = 0.1. Our curve is similar to a scaling 
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curve of Y vs X which is predicted by Takayama and 
Hukushima 8 from the Monte Carlo simulation on the 3D 
Ising EA model under the iTshift protocol. The scaling 
curve is described by 

Y = X-cX 1+1 / 5 , (16) 

for X <C 1, where c is a constant. In Fig. ITTT b). for 
comparison we show the plot of Y vs X given by Eq. ([THjl 
with 5 = 0.81, when c is changed as a parameter. Our 
data are fitted well with Eq. (|T?)J) with c = 0.15 for X < 
0.3. Note that when c and S are adjustable parameters, 
the least squares fit of the data (Y vs X) (0 < X < 0.6) 
to Eq. (|TB|) yields the parameters c = 0.71 ± 0.06 and 
S = 0.81 ±0.07. 

In summary , we find a scaling relation that all the set 
of Y = R cr /RH vs X = Rw/Rh obtained under the H- 
shift protocol fall well on a single scaling curve. The curve 
of Y vs X clearly shows a gradual dynamic crossover from 
the relation Y = X for X <C 1 to Y « 1 for large X. In 
this sense, the iJ-shift aging process is nothing but a 
dynamic crossover from the SG state in H = to the 
paramagnetic state in H > in the equilibrium limit. 

What is the relation between H and T, depending on 
the value of X c ? Note that X c ^0.1 for the deviation of 
Y from the linear relation Y = X, and X c « 1.6 for the 
saturation of Y. From the relation X = X c in Eo. (|15a|) . 
the magnetic field H can be derived as 

H « XyV /5 (l - T/T c ) a ^f(t w /rS)- {b/5)T/Tc . (17) 

In Fig. H2T a) we show the plot of H vs T defined by 
Ea.(|17|) as X c is varied as a parameter (0.5 < X c < 
1.5), where t w = 1.0 x 10 4 sec, r * = 5.27 x 1(T 6 sec, 
C = 1.5 x 1(T 4 , S = 0.81, b/S = 0.198, and a eff = 
0.55. It seems that the H-T diagram for the fixed X c 
is similar to Eq.fTJ with p = 3/2 (the AT prediction). 
In fact, the least squares fit of these data of H vs T 
to Eq.© with t w = 1.0 x 10 4 sec for 2 < T < 3.6 K 
yields the parameters A and p. The value of p is equal 
to 1.50 ± 0.03, independent of X c , while the value of A 
increases as increasing X c : A = (8.20 ± 0.15)X 1 - 31±0 06 
[kOe]. 

Finally we discuss the nature of the line T g o(H) 
(dAx/dT vs T) shown in Fig. [3 This line is well de- 
scribed by Eq.© withp = 1.52±0.10and A = 3.38±0.67 
kOe. If the above model is applicable to this line, 
the value of X c is estimated as X c = 0.50. Figure 
H2T b) shows the plot of H vs T, where H is calculated 
from Ea.([T7|) for each T. Here we choose tw = 100 
sec which is appropriate for the measurements of Xzfc 
and Xfcs and X c = 0.28 which is rather smaller than 
X c = 0.50. We find that the line T g0 (H) (dAx/dT vs 
T) for 2.5 < T < 3.6 K lie well on the prediction from 
Ea.(|17j). This result also suggests that the line T g o(H) 
(dAx/dT vs T) is also explained in terms of the droplet 
picture. 
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FIG. 12: (Color online) (a) H-T line (denoted by the solid 
lines). The value of H is calculated from Ea. (fTTJ) with tw = 
1.0 x 10 4 sec and X c being varied as a parameter (0.5 < X c < 
1.5), where r * = 5.27 x 10~ 6 sec, ( = 1.5 x 10~ 4 , 5 = 0.81, 
b/5 — 0.198, and a e // = 0.55. The dotted lines denote the 
least-squares fitting curves of the H vs T relation for each X c 
to Eq.JI} with p w 1.5. (b) H-T line. The value of H as a 
function of T is calculated from Ea. (|17f> with X c = 0.28 and 
tw — 100 sec, where the values of Tq , £, 5, b/5, and a e ff are 
the same as those for (a). The line T g o(H) (dAx/dT vs T) 
(•) of Fig.|5]is also shown for comparison. 



VI. CONCLUSION 

Cuo.sCoo.sC^-FeCls graphite intercalation compound 
is a three-dimensional short-range Ising SG with a spin 
freezing temperature T c (= 3.92 ±0.11 K). The stability 
of the SG phase in the presence of a magnetic field H has 
been studied from the spin freezing temperature Tf(oj, H) 
and the peak time t cr (T, H) of S{t) vs t for T < T c . We 
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find that both Tf(uj, H) and t cr (T, H) show certain scal- 
ing behavior predicted from the droplet picture. The 
ratio Y (= Rcv/Rh) is well described by a scaling func- 
tion of X (= Rw/Rh) in the form of Y = X - cX 1+1 / 5 
with c = 0.71 ±0.06 and S = 0.81 ±0.07. This result indi- 
cates the instability of the SG phase in the thermal equi- 
librium in a finite magnetic field: a dynamic crossover 
occurs from SG dynamics to paramagnetic behavior. 
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